Direct transition from quantum escape to phase diffusion regime in YBaCuO 

biepitaxial Josephson Junctions 
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Dissipation encodes interaction of a quantum system with the environment and regulates the 
activation regimes of a Brownian particle. We have engineered grain boundary biepitaxial YBaCuO 
junctions to drive a direct transition from quantum activated running state to phase diffusion regime. 
The cross-over to the quantum regime is tuned by the magnetic field and dissipation is encoded in 
a fully consistent set of junction parameters. To unravel phase dynamics in moderately damped 
systems is of general interest for advances in the comprehension of retrapping phenomena and in 
view of quantum hybrid technology. 

PACS numbers: 05.40.Jc, 74.50. +r, 85.25.Cp 
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The impressive development of superconductive sys- 
tems in the field of quantum information processing and 
the expertise gained on manipulating coherent entangled 
states and different coupling regimes with the environ- 
ment *j] have boosted research on several complemen- 
tary aspects of coherence and dissipation. Due to their 
design scalability and to the flexibility in controlling the 
level of damping, Josephson systems have proven to be 
a fantastic test-bench for studying fundamental physics 
problems such as the quantum superposition of alive 
and dead states of a Schrodinger's cat[2(, the behavior 
of an artificial atom in cavity quantum electrodynamics 
experiments Q, or measurements of quantum coherence 
in macroscopic systems 0, HI- 

Both the time evolution of the position of a Brownian 
particle and the electrodynamics of a Josephson junction 
(J J) [1,0 can be described by a Langevin equation [8] of 
the form 



<p + tp/Q + dU/d<p = £(t) 



(1) 



In this equation the time is normalized to l/u p , with 
Up = (2eI co /hC) 1 ^ 2 representing the plasma frequency 
at zero bias current and I co and C being the junction 
critical current in absence of thermal fluctuations and 
capacitance, respectively. The potential U is the well 
known periodic "washboard" potential associated to the 
dynamic of a JJ, U(ip) = —Ej^co&ip+Y-ipj, where 

tp(t) is the superconductive phase and Ej = hI co /2e is 
the Josephson energy. £(t) is a white noise driving force 
such that 



<€(*)) =0; <£(*),£(*')> = y/k B T/QEj6(t-f). (2) 

The parameter Q is given by Q = u p RC 0, 0] with R 
the shunt resistance. In a more general approach Q has a 
frequency dependence [§] better responding to the need of 



including external shunting impedance. In our case the 
simplest approach based on the Q factor calculated only 
at plasma frequency (loi gives a satisfying account for 
experimental data, as demonstrated below. It is a natural 
choice to use a JJ in order to study the Brownian motion 
of a particle in a dissipative tilted periodic potential (see 
FigQH,b). 

In this letter we demonstrate a direct transition from a 
running state, obtained following a quantum activation, 
to diffusive Brownian motion in YBaCuO JJs. Multiple 
retrapping processes in subsequent potential wells char- 
acterize phase re gim es where diffusive phenomena play 
a relevant role [9l4l3j|(see FigQJ)). The relevant parame- 
ters driving the occurrence of these phenomena are the 
operational temperature T, the damping factor Q and 
the critical current I co . The various operation scenar- 
ios for a JJ have been condensed in a phase diagram 
by Kivioja et al. fioj] who have shown that by spanning 
the (EjjksT) parameter space it is possible to engi- 
neer all different regimes ranging from phase diffusion 
and thermal activation to macroscopic quantum tunnel- 
ing (MQT). MQT takes place not only for low values of 
dissipation (Q >> 1), but also for intermediate levels of 
dissipation (1 < Q < 5). We explore a new region of this 
phase diagram, made available by the different ranges of 
I co and of the standard deviation of the switching dis- 
tribution a offered by these junctions when compared 
with most low temperature superconductors (LTS) JJs. 
The moderately damped systems are particularly signifi- 
cant and promising to address and quantify interactions 
of a quantum system with the environment which 
is, apart from the its intrinsic interest for fundamental 
physics, a cornerstone for the development of whatever 
quantum hybrid technology. 

The experimental observation of the various regimes 
in a JJ is based on the measurement of the switching 
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FIG. 1. a) Thermal (red dashed line) or quantum activated escape in the tilted periodic potential. Quantum escape is 
represented both for very low ideal (Q >> 1 - continuous black line) and high (1 < Q < 5 - dashed black line) levels of 
dissipation respectively, b) Diffusive motion due to multiple escapes and retrapping in subsequent potential wells, c) Schematic 
escape rates of MQT, TA and retrapping processes at T < T x (top panel) and T > T x (bottom panel) d) Sketch of the off-axis 
biepitaxial junctions. The JJ is formed at the boundary between (001) YBaCuO and (103) YBaCuO electrodes. Three interface 
orientations 9 — 0°, 60° and 90° are shown as examples. Different interface orientations (lobe vs. lobe, node vs. lobe and any 
configuration in between) can be achieved with proper patterning of the seed layer [23j| . e)Picture of the device along with a 
block diagram of the experimental setup. In the bottom part of the image (in blue) is the YBaCuO (001) electrode while in 
the top part (gray) the needle-like YBaCuO (103) grains are visible. 



current distribution (SCD) and the study of the behavior 
of its first and second momenta (the mean I and the 
width er) as function of temperature. In an underdamped 
junction (Q > 10) [la ], below a crossover temperature T x 
the escape process is mostly due to MQT (see Fig. [I]}), 
marked by a temperature independent a, while above T x 
the process of escape is due to thermal activation (TA) 
above the potential barrier, with a distinctive increase 
of a with temperature. In moderately damped junctions 
[l-Qil [H Illm8t with 2 < Q < 5 a transition from TA 
to 'phase diffusion' (PD) regime occurs at a crossover 
temperature T* > T x . T* marks a distinctive change in 
the sign of the temperature derivative of a, with da/dT > 

for T < T* and da/dT < for T > T* . When escape 
out of a well occurs at too low currents in PD regime [9J, 
the energy gained by passing from one well to the next 
one barely exceeds the dissipative losses and the particle 
eventually gets retrapped, basically diffusing to the very 
next wells. The SCD histograms move to lower currents 

1 till they touch the limit In = (A*I co /-k)\/Q in the most 
common cases with Q ;§> 1 [3, [l9| . 

We report measurements of SCD in a temperature 
range from 20mK to 2.2K. Our data are characterized 
by two distinct regimes. Below 135mK the widths of 
the SCDs show no significative variation. This is a typ- 
ical signature of a quantum activation regime. Above 
135mK the negative temperature derivative of a is con- 
sistent with a diffusive motion due to multiple retrapping 



in the potential wells. This regime has been fitted using 
existing theories on phase diffusion (l8j with a damping 
factor Q = 1.3. The observations are qualitatively con- 
sistent with the escape rates sketched in FigJTb. We have 
substantially engineered a device with T* <T X . For tem- 
peratures T well below T x , MQT contributions to escape 
rates are larger than those coming from both thermal es- 
cape and multiple retrapping processes (Fig(TJ: top panel) 
differently from the case T > T x (FigJTJ; bottom panel). 
The contiguity between quantum escape (T < T x ) and 
phase diffusion (T > T x ) leads to MQT phenomena char- 
acterized by low Q values and not necessarily to quantum 
phase diffusion. This MQT process can be represented in 
Fig. |TJi as a dashed line to manifest interaction with the 
environment [2lj |. and responds to what it can be possi- 
bly experimentally ascertained. This phenomen ology i s 
quite distinct from all previous studies 

EE lid-lid . 

where in the transition to quantum activation, retrapping 
processes decay faster than thermal escape, and from the 
work of Yu et al. (20l | , where the occurrence of a quantum 
activated phase diffusion has been claimed. In ref [20j j . 
the semiclassical nature of their quantum phase diffusion 
is testified by the dependence of a on the temperature 
over the entire temperature range, and the transition is 
as a matter of fact revealed by a change of the temper- 
ature derivative of a [2(j. MQT processes are substan- 
tially followed and assisted by thermally-ruled retrapping 
processes. However, a fully quantum account of phase 
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FIG. 2. a) Measured switching current probability distribution P(I) at different bath temperatures, for sample A. The inset 
shows the device current voltage characteristic measured at 30mK. The reference value for the threshold detector is also 
displayed, b) Temperature dependence of the standard deviation, a, of the switching distributions for sample A. The dash- 
dotted line marks the temperature-independent SCD widths in the quantum tunneling regime, the red solid line is the results 
of simulations in the diffusive regime with a damping parameter of Q=1.3. The inset shows temperature dependent data for 
sample B acquired at two different values of the applied magnetic field, c) Simulated thermal behavior of the width of the 
switching histogram for several values of the Q damping parameter. In the inset we report the dependence of the turn-over 
temperature T* on the damping parameter. 



fluctuations passes through the 'empirical' condition of 
a Josephson energy much larger than Coulomb energy, 
Ej » E c with Ec = e 2 /2C (see below), given by Ian- 
siti et al. |22j. This condition does not occur both in 
this experiment and in the work of Yu et al. |20| , which 
both represent complementary significant advances to- 
wards the observation of a fully-quantum phase diffusion, 
better defining its domain through the border with com- 
peting processes. 

We have used YBaCuO off-axis GB biepitaxial JJs, 
whose scheme is shown in Fig. [TJi [231425 1. The 
GB is determined at the boundary between the 103- 
oriented grains growing on the bare substrate and the 
001 grain growing on the CeO^ seed layer. We have 
engineered junctions on (Lao. a^o. 7) (^-^0.65^00.35)03 
(LSAT) rather than SrTi0 3 (STO) substrates, where 
MQT in a high temperature superconductor JJ was first 
demonstrated (SB)]. The new design fully responds to the 
need of reducing stra y ca pacitances and of better isolat- 
ing the GB behavior [26[. Specific capacitances are one 
order of magnitude lower than those measured on STO- 
based devices [2^, [26J. Dynamical junction parameters 
can be tuned by choosing the interface orientation indi- 
cated by the angle 8 in Fig. [TJi, which also sets d-wave 
induced effects (23, HI- 

To study the escape rates of YBaCuO Josephson junc- 
tions we have thermally anchored the sample to the mix- 
ing chamber of a He3/He4 Oxford dilution refrigerator 
and performed measurements of the junction switching 
current probability. A schematic representation of our 
experimental set-up is shown in FigJT^. A full descrip- 
tion of the apparatus is discussed in detail elsewhere [l7| . 
Filtering is guaranteed by a room temperature electro- 
magnetic interference filter stage followed by low pass RC 



filters with a cut-off frequency of 1.6MHz anchored at 1.5 
K, and by a combination of copper powder and twisted 
pair filters thermally anchored at the mixing chamber of 
the dilution refrigerator. The bias current of the junction 
is ramped at a constant sweep rate dl/dt — 17.5/iA/s 
and at least 10 4 switching events have been recorded us- 
ing a standard technique. These measurements, collected 
over a wide range of temperatures are reported in FigJ^i. 
A progressive broadening of the histograms occurs when 
lowering the temperature, which is a distinctive feature 
of the phase diffusion regime For tem- 

peratures below about 135mK, the histograms overlap. 

The temperature dependence of the width a of the 
SCD curves is shown in Fig. As a test of fidelity for 
our fabrication process and our experimental setup we 
report data for two different samples with interface ori- 
entations of 75° for sample A and 50° for sample B. The 
junctions are 2.5 /j,m and 2.0 /im wide respectively, and 
the film thicknesses are 100 nm and 250 nm. We have 
selected interface orientations with robust overlap of the 
d-wave lobes on both sides of the junctions[24| and Ej 
of the order of 3 meV. For both samples above a tem- 
perature T* the data display a decrease in the width of 
the distribution in agreement with a process of multiple 
escapes and retrapping typical of the diffusive motion. 
The high temperature region (135mA' < T < 2K) allows 
a reliable estimation of the damping parameter. Simula- 
tions for Q = 1.30 are reported as a red line in Fig. (2b 
resulting from the integration of the Langevin equation 
([T]) with a Bulirsh-Stoer integrator using a noise affected 
Runge-Kutta Monte Carlo algorithm [17j. 

In Fig. [2k simulated thermal behavior of a is reported 
for different values of the Q damping parameter ranging 
from 1.2 to 5. For each of these curves T* approximately 
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TABLE I. Device parameters. 



Sample 


H (G) 


I CO (M) 


R (fi) 


C (fF) 


Q 


T x (mK) 


A 





1.20 


84 


64 


1.30 


135 


B 





1.79 


64 


74 


1.28 


144 


B 


12 


1.42 


64 


74 


1.14 


122 



indicates the transition temperature from thermal activa- 
tion to the diffusive regime. Q tunes T* as shown in the 
inset of Fig. [2fc and modifies the slope of the <x(T) fall-off 
at higher temperatures. The capability to numerically re- 
produce this region makes it possible to estimate Q with 
high precision. In our case Q = 1.30±0.05 closely fits the 
data and determines a T* value not larger than 100 mK. 
The section below T* faithfully reproduces the expected 
T 2 / 3 dependence for a thermally activated regime 
(solid line) as an additional test of consistency. In Fig. 

the MQT section is missing. It would attach below T x 
to each of the curves with its characteristic saturation in 
tr, as shown in Fig. [2Jd in fitting experimental data. 

In facts below 135 mK a is almost independent of tem- 
perature and consistent with an escape dominated by 
quantum tunneling. Thus 135 ± 10 mK represents the 
cross-over temperature T x providing, once known the Q 
value from the phase diffusion fitting, first C = 64 fF 
and the plasma frequency uj p ~ 38 GHz. This value of 
the capacitance C is also in agreement with the value 
extrapolated using the estimated specific capacitance 
C s ~ 2 • 10" 5 Fcrnr 2 for YBaCuO BP JJs on LSAT sub- 
strates 24 , [26} . As a consequence Ej ~ ZmeV results 
to be much larger than E c ^ SfieV (Ej/E c ^ 1000) 



and the system is in conditions nominally far from those 
which are considered to be favorable for the observation 
of quantum phase diffusion [22] . 

In analogy to what commonly done to prove MQT in 
underdamped i unctions [lB|. we use the magnetic field to 
tune in situ the junction parameters and T x to unambigu- 
ously prove MQT as source of the saturation of a below 
T x . In the inset of Fig. [2p we report the temperature 
dependence of a measured for sample B at two different 
magnetic fields of 0G and 12G respectively. H=12G low- 
ers the critical current I co reducing at the same time the 
quantum crossover temperature. Relevant device param- 
eters are summarized in Table [U 

In Fig. [3] we report a (Q, ksT/Ej) phase diagram, 
which summarizes the various activation regimes (Ioj|. 
The transition curve between the phase diffusion regime 
and the runnin g st ate following thermal (experimentally 
observed in I(J 13, [H, or quantum (experimentally 
observed in this work) activation has been determined 
numerically by varying the damping factor Q as func- 
tion of the ratio between the thermal energy and the 
Josephson energy. The filled circles are experimental 



data obtained by fitting the escape rates T as a func- 
tion of the ratio between the barrier height and the es- 
cape energy, u = AU/kBT esc [28[ (shown in the inset of 
Fig. [3]). The obtained values fall within the region of 
the diagram that displays a direct transition from phase 
diffusion to quantum activation. The escape rates have 
been calculated from the switching distributions using a 
standard procedure [29]. In the quantum activated regime 
the switching distributions are asymmetric and skewed to 
the left, and T values all fall onto the same line, as it is 
the case for the reported data from T=30mK to 108mK. 
Retrapping processes cause a progressive symmetrization 
of the switching distribution which translates into a de- 
viation of the experimental escape rates from the ideal 
exponential behavior 17 , Tj^ . 
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FIG. 3. (QyksT I Ej) parameter space, showing the various 
activation regimes. The transition curve between the PD 
regime and the running state has been extrapolated through 
numerical simulations, the sideband curves mark the uncer- 
tainty in our calculation and are due to the temperature step 
size. The crossover temperature between the various regimes 
scales with Ej. The data points refer to the Q values of our 
sample and show the direct transition from quantum acti- 
vated running state to the PD regime. In the inset we show 
the experimental escape rates (symbols) as function of barrier 
height to escape energy ratio along with the theoretical fits at 
different T. The Q values used for the fits are the same shown 
in the phase diagram. 



In summary, by exploring a new region of the 
(Q, fcsT/ Ej) phase diagram we have demonstrated a 
direct transition from quantum activation to diffusive 
Brownian motion in GB YBaCuO JJs. This experiment 
sets another milestone in the study of the influence of dis- 
sipation on the switching statistics of JJs 12 , T^, 21 , 30l — 
|32j demonstrating novel balancing between MQT, ther- 
mally activated and retrapping escapes, thus paving the 
way to the observation of fully-quantum phase diffusion, 



5 



and is of particular relevance to understand interaction of 
a quantum system with the environment [14|. The com- 
bined experimental and numerical investigation here pre- 
sented has the potential to offer a new tool to study the 
interplay between coherence phenomena and dissipation 
down to the quantum regime in a wide variety of systems. 
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